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We analytically compute the spectrum of the spin zero glueballs in the planar limit of pure 
Yang-Mills theory in 2-1-1 dimensions. The new ingredient is provided by our computation of a new 
non-trivial form of the ground state wave-functional. The mass spectrum of the theory is determined 
by the zeroes of Bessel functions, and the agreement with large N lattice data is excellent. 


The understanding of the non-perturbative dynamics 
of Yang-Mills theory is one of the grand problems of the¬ 
oretical physics. In this letter we announce new analyt¬ 
ical results pertaining to the spectrum of the spin zero 
glueballs of 2 -I-1 dimensional Yang-Mills theory [ij . This 
theory is expected on many grounds to share the essential 
features of its 3 -I- 1 dimensional cousin, such as asymp¬ 
totic freedom and confinement, yet is distinguished by 
the existence of a dimensionful coupling constant. Here, 
we determine the ground state wave-functional in the pla¬ 
nar limit, the knowledge of which enables us to do the 
necessary computations regarding the mass gap, string 
tension and the glueball spectrum. The results are in 
excellent agreement with the lattice data in the planar 
limit 1^. Full technical details will appear in a longer 
publication Q. 

Our approach, as in our previous work 0, is based 
on the remarkable work of Karabali and Nair [^. The 
Karabali-Nair approach can be summarized as follows. 
Consider an SU{N) YM 2+1 in the Hamiltonian gauge 
Aq = 0. Write the gauge potentials as Ai = 
for i = 1,2, where are the Hermitian N x N matri¬ 
ces in the SU{N) Lie algebra with the 

normalization 2Tr{t°‘t^) = 5“^. Define complex coordi¬ 
nates z = xi — ix 2 and z = Xi + 1 x 2 , and furthermore 

2A“ = A? -b iA^, 2A“ = A? - iA^. 

The Karabali-Nair parameterization is 

A = -d,MM-\ A =(1) 

where M is a general element of SL{N, C). Note that a 
(time independent) gauge transformation A —> gAg~^ — 
dg A gAg~^ - dg where g S SU{N), be¬ 
comes simply M —> gM. Correspondingly, a local gauge 
invariant variable is H = M. The standard Wilson 
loop operator may be written 

<i)(C) = TrPexp{j) dz 


The definition of M implies a holomorphic invariance 
M{z,z) M{z,z)h^{z) 

(3) 

Mt(z,z) ^ h{z)M\z,z) 

where h(z) is an arbitrary unimodular complex matrix 
whose matrix elements are independent of z. This is dis¬ 
tinct from the original gauge transformation, since it acts 
as right multiplication rather than left and is holomor¬ 
phic. Under the holomorphic transformation, the gauge 
invariant variable H transforms homogeneously 

H{z,z)^h{z)H{z,z)h\z). (4) 

The theory written in terms of the gauge invariant H 
fields will have its own local (holomorphic) invariance. 
The gauge fields, and the Wilson loop variables, know 
nothing about this extra invariance. We will deal with 
this by requiring that the physical state wave functionals 
be holomorphically invariant. 

One of the most extraordinary properties of this pa¬ 
rameterization is that the Jacobian relating the measures 
on the space of connections C and on the space of gauge 
invariant variables H can be explicitly computed [Sj 

dpL[C] = (5) 

where ca is the quadratic Casimir in the adjoint represen¬ 
tation of SU{N) {cA = N), cr is a constant determinant 
factor and 

Swzw{H) = -^J d^z TrH-^dHH-^dH 

+-^ J d^x e>^^^TrH-'^d^,HH-^dyHH-^dxH 

is the level— CA SU{N) Wess-Zumino-Witten action, 
which is both gauge and holomorphic invariant. Thus 
the inner product may be written as an overlap integral 
of gauge invariant wave functionals with non-trivial mea¬ 
sure 

( 1 | 2 ) = / dp,[H]e^‘^^^'^^'^^^'>AllAi2. 


( 2 ) 


(7) 


2 


We note that in this norm, 'I' = 1 is normalizable. 

From these expressions it is clear that a useful gauge- 
invariant variable is the current 

J = 2Ad,HH-^ (8) 

TT 

which transforms as a holomorphic connection 


The standard YM 2+1 Hamiltonian 


T~(-ym — I Tr [ "I— 2 — 

9ym 


( 10 ) 


J ^ hJh ^ + —dzhh (9) 

TT 


Note that dJ transforms homogeneously, and a 
holomorphic-covariant derivative is given by 


CA 


can be also explicitly rewritten in terms of gauge invari¬ 
ant variables. The collective field form of this Hamil¬ 
tonian (which we will refer to as the Karabali-Nair Hamil¬ 
tonian) can be easily appreciated from its explicit form 
in terms of the currents as follows 


nKN[J] = T + v = 


rix) 


SJ°-{x) 


y) 


lx,y 


6J°-{x) SJ^(y) J mcA 


Brar 


( 11 ) 


where 


_ 9 yM^A y. _ CA Sab _ i fabcJ^^jx) 

2tT ' ab , ^2 _ y ^2 ^ — y) 

( 12 ) 

Interpreted as a collective field theory, one can expect 
to compute, at large N, correlators of gauge invariant 
operators. Note that the magnetic field is 

27r - 

B = -2 M'<-^d{dHH-^)M'< = -(13) 

CA 


The derivation of this Hamiltonian involves carefully reg¬ 
ulating certain divergent expressions in a gauge invariant 
manner 0, 0- We note that the scale m is essentially 
the’t Hooft coupling. 

The purpose of this letter is to determine masses of 
some of the lowest lying glueball states. To do so, we wish 
to determine the form of the vacuum wave functional and 
make use of the planar limit. Accordingly, we take the 
following ansatz for the vacuum wave functional 


This form of the wavefunctional is explicitly gauge and 
holomorphic invariant. The kernel AT is a formal Tay¬ 
lor expansion of L = {DB + BD)/2m^, while the ellipsis 
contains terms higher order in dJ (or B). This wave- 
functional has the form of a “generalized coherent state” 
appropriate to large N 0 , but its form is not completely 
dictated by large N counting. The form of the ansatz, as 
we shall see, is sufficient to capture the mass spectrum 
of gauge invariant states, which we will probe using local 
operators. The large N limit ensures that these states 


are non-interacting, but we are also neglecting the size of 
the states by using local probes. (For further details on 
these points, see Ref. 0 .) 

In order to be physically sensible, K should have cer¬ 
tain properties at long and short distances. We derive 
these properties below. In particular, the low momen¬ 
tum (large ’t Hooft coupling) limit, m?, of the 

vacuum wave functional is easily determined to be of the 
form 

(Equivalently, at low momentum, we should have K —> 
1.) This wavefunctional provides a probability measure 
'Fq^'o equivalent to the partition function of the Eu¬ 
clidean two-dimensional Yang-Mills theory with an effec¬ 
tive Yang-Mills coupling = ^9ym- Using the results 
from 0, Karabali, Kim and Nair deduced the area law 
for the expectation value of the Wilson loop operator 


($) ^ exp(— ctA) (16) 

with the string tension following from the results of 0 


4 N^-1 

<^-9ym—^ 


(17) 


This formula agrees nicely with extensive lattice simula¬ 
tions 0 , and is consistent with the appearance of a mass 
gap as well as the large N ’t Hooft scaling. 

Coming back to the derivation of the vacuum wave 
functional, we argue in Ref. 0 that operators On = 
f BJL'^BJ, which would appear in a series expansion of 
K{L), satisfy 


TOn — (2 -l- n)mOn -I- ... 


(18) 
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In Ref. [3] (see Sec. 3 and App. A), we have presented 
a series of calculations supporting this important result. 
Further evidence is provided by lattice considerations. 
Given this, we can formally write 

TKiL)-^^-^[L^K{L)]. (19) 

The full vacuum Schrodinger equation, combining all 
contributions self-consistently to quadratic order in dJ 


'Hkn'^o = Eq'^’o = 


trdJTZdJ ' 


^ 0 , ( 20 ) 


with suitable subtractions, then formally leads to the fol¬ 
lowing differential equation for K 

^7^ = -iF-^^[i^(L)]+LiF" + l = 0. (21) 

TT 2 clL 

In this equation, the final term is the contribution of 
the potential term of the KN Hamiltonian, while the 
second to last term arises from the fl-term in the kinetic 
energy. Eq. m comes by consistently keeping all terms 
quadratic in dJ in the Schrodinger equation. 

Although this equation is non-linear, it is easily solved 
by substituting K = —y'j2y; the resulting equation may 
be recast as a Bessel equation. The only normalizable 
solution has the correct physical asymptotics for large and 
small L and is given by 


x(L)= 1 

/LJi(4VI) 


( 22 ) 


where Jn denotes the Bessel function of the first kind. 
This remarkable formula encodes information on the 
spectrum of the theory, as we show below. We note 
that this kernel has the following asymptotics (where 
L ^ —p’^l^m?) 


p ^ 0, K ^ 1] p^oo, K ^ 2m/p (23) 


consistent with confinement and asymptotic freedom, re¬ 
spectively. 

In order to determine the spectrum, we factorize suit¬ 
able correlation functions at large distances. The op¬ 
erators appearing in the correlation functions will have 
definite quantum numbers, which will be inherited 
by the single particle poles contributing to the correlation 
function. 

As a first example, we consider the states which 
may be probed by the operator TrdJdJ. We have 

{TridJdJ)^ TridJdJ);, ~ {K-\\x - y\)f . (24) 

Here, we have computed the correlation function in the 
planar limit given our knowledge of the vacuum wave- 
functional. 


TABLE I: glueball masses in YM 2 + 1 . All masses are in 

units of y/a. AdS/CFT computations!^ are also given for 
comparison. The percent difference between our prediction 
and lattice data is given in the last column. 


State 

Lattice, N —> 00 

Sugra 

Our prediction Diff, % 

0 ++ 

4.065 ± 0.055 

4.07(mput) 

4.10 

0.8 

0 ++* 

6.18 ±0.13 

7.02 

5.41 

12.5 

0 ++** 

7.99 ± 0.22 

9.92 

6.72 

16 

Q+ + .** 

9.44 ± 0.38“ 

12.80 

7.99 

15 


“Mass of 0++*** state was computed on the lattice for SU{2) 
only i. The number quoted here was obtained by a simple rescal¬ 
ing of SU{2) result. 


To proceed further, we note the identity 
_ 2 :/ \" 1 


(25) 


where jv,n are ordered zeros of the Bessel functions. For 
example, the first few zeros 0 of J 2 {z) are ^ 2,1 = 5.14, 
32,2 = 8.42, j 2,3 = 11.62, , 72,4 = 14.80, etc. Apart from 
additive constants, we then deduce 


K-^{k) 


1 

2 


E 


M^ + k^ 


(26) 


where Mn = j 2 ,nt 7 i/ 2 . The Fourier transform at large 
|a; - y\ is 


K {\x-y\) = - 


— , y (M„)3/2e-“"l^-!^l. 

4^2TT\x-y\ 

(27) 

In particular the correlator mentioned above is 


1 

327 r|a; 


y (28) 

n, m—1 


Note that each term here has the correct \x — y\ depen¬ 
dence for a single particle pole of mass -|- Mm in 2 -|-1 
dimensions. The 0'^“*' glueball masses are: 


Mq+ + 

= Ml 

+ 

Ml = 

= 5.14m 

Mo++» 

= Ml 

+ 

M 2 = 

= 6.78m 

Afo++»* 

= M 2 

+ 

M 2 = 

= 8.42m 

Mg + + ../ 

= Ml 

+ 

M 3 = 

= 8.38m 

Mo++*** 

= M 2 

+ 

M 3 = 

= 10 . 02 m, 


Since m is not a physical scale, we should re-write these 
results in terms of the string tension. Given equations 
presented above, at large N we have ~ Our 

results are given in Table 1. Several comments are now 
in order. First, note that we have been able to predict 
masses of the 0 +“'' resonances, as well as the lowest lying 
member, in contrast to the original results of Karabali 
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TABLE II: 0 glueball masses in YM 2 + 1 . Columns are as 
in Table I. 


State 

Lattice, N ^ 00 

Sugra 

Our prediction 

Difl,% 

0-- 

5.91 ±0.25 

6.10 

6.15 

4 

0““* 

7.63 ± 0.37 

9.34 

7.46 

2.3 

0 

1 

1 

8.96 ± 0.65 

12.37 

8.77 

2.2 


and Nair (which differ significantly numerically). The 
supergravity results listed in the table are a result of cal¬ 
culations 0 using the AdS/CFT correspondence 0 ; in 
that case, the overall normalization was not predicted but 
was determined by fitting to the lattice data, for example, 
to the mass of the lowest O'*'’*' glueball. Our results for 
the excited state masses differ at the 10-15% level from 
the lattice results. We note that precisely these masses 
are more difficult to compute on the lattice , and thus 
the apparent 10 — 15% discrepancy may be illusory. [l^ 
Finally, we note that there is some interesting approxi¬ 
mate degeneracies in the spectrum. 

Let us move on to a discussion of the 0 glueball 
resonances. In this case, our predicted masses are much 
closer to the lattice data, which we believe to be more 
reliable in this case. We may probe these states with 
the operator Tr dJdJdJ . We are thus interested in the 
correlation function m 

{TiidJdJdJ)^ Tr{dJdJdJ)y) ~ {K-^{\x-y\)f . 

(30) 

Using the results given above, we obtain glueball masses 
which are the sum of three M„’s. 

Mq-- = Ml-b Ml-b Ml = 7.70 to 

Mo—. = Ml-b Ml-b M 2 = 9.34to (31) 

Mq--.. = Ml -b M 2 -b M 2 = 10.99to. 

These results are compared to lattice and supergravity 
data in Table II. We see that the resulting masses are 
within a few percent of the lattice data, and are much 
better than the supergravity predictions. 

Our results suggest that there exist hidden constituent 
as well as integrable structures in 2-b 1-dimensional Yang- 
Mills theory. Note that the full integrability of 2 -b 1- 
dimensional pure Yang-Mills has been suspected for some 
time 0 In a longer paper we will more carefully 
explain our techniques and results and will investigate 
other glueball states and the corresponding Regge 
trajectories. It has not escaped our attention that a sim¬ 
ilar parameterization may be used in 3-bl dimensions in 
a variational sense and preliminary numerical results are 
encouraging. 
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